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Abstract: 

Many experimental evidences indicate the presence of a ionizing background radiation flux at large 
redshifts, whose nature is doubtful. A lot of informations about the characteristics of such a background can 
be obtained both from the study of the Gunn Peterson effect and from the so-called "proximity effect" . 

In some previous works I suggested the possibility that this ionizing flux comes from the quantum 
evaporation of primordial black holes (PBHs): here, I discuss the constraints that the experimental mea- 
surements put upon the free parameters of this reionization model and I try to verify its reliability. In 
particular, the radiation intensity of the background at the hydrogen Lyman edge, as inferred from the 
proximity effect, enables me to determine an upper limit to the PBHs average relics mass; due to our poor 
knowledge of the ultimate fate of the evaporating black holes, this limit represents an important theoretical 
information. In the second part of this paper I study the absorption of the ionizing background due to Lya 
clouds: in particular, I discuss this phenomenon in presence of different absorption levels and I calculate the 
HI Gunn Peterson optical depth tqp{z); from a comparison with the experimental data of Giallongo et al. 
{'^GP. HI < 0.02 ± 0.03) I obtain a constraint on the intergalactic medium density parameter, namely ^igm 
< 0.020. 

A study of the characteristics of the absorbers is also performed: I determine the hydrogen gas density 
nH,c and the column density Nhi for Lya clouds; a satisfactory agreement with the available experimental 
data is obtained in the case of expanding, adiabatically cooled clouds. Finally, the same kind of analysis 
is performed for He II: in this case, the theoretical optical depth I obtain is smaller than the preliminary 
experimental lower limit of Jakobsen et al. {tgp > 1.7). 
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1. INTRODUCTION 

The Gunn Peterson effect [1] is known since 1965 as a test of our knowledge of the ionization history 
of the Universe: the presence of an uniform distribution of neutral hydrogen at a redshift z should indeed 
be proved by the presence of an absorption trough in the bluesidc of the Ly^ emission spectrum of the 
observed quasars; in fact, the Lya resonance absorption line, (A=1216 A), is shifted by a factor (1 + z) 
into the visible, blue part of the spectrum. However, no trough is experimentally observed, in spite of the 
presumable sharpness of this effect; this fact probably reveals our unsatisfactory knowledge of the status of 
the Universe during the post-recombination era. 

On the basis of the pioneering Schmidt's observations [2] of the quasar 3C9 at 2; = 2.016, Gunn and 
Peterson [1] obtained an integrated Ly^ optical depth tgp < 0.5, corresponding to an intergalactic neutral 
hydrogen number density 

UHiiz) = rGp{z){l + z){l + 2q,zfl' (^^^^^) = 

^ TGpjz) (l + z) {l + 2qozy/\ 
4.14 X IQio /i-i 

In cq. (1.1), Ho is the present value of the Hubble constant, equal to 100 h km s^^ Mpcr^, 0.4 < h < 1, 
qo = 1/2 is the deceleration parameter, Ua = 2.46 x 10^^ sec^^ is the HI Lya frequency, rup the proton mass 
and / = 0.416 is the oscillator strength corresponding to the Lya transition [1]. From eq. (1.1), nHi at 
z = 2 is equal to 6 x 10""'^'^ h cm~'^; subsequent works [3] still lowered this limit, that now approximately 
is 15 times smaller than the Gunn Peterson original estimate; for instance, Steidel and Sargent found [4] 
nHi{z = 0.0) < 9.0 X 10-1^ h cm-^. 

These values are smaller that the cosmic abundances: the inobservance of the Ly„ absorption trough is 
probably attributable to a high ionization level of the Universe rather than to an effective paucity of neutral 
hydrogen. 

Many causes of such a ionization have been postulated: for instance, shock heating phenomena [5], [6], 
[7] , high mass stars in primordial galaxies and unseen quasars, hidden due to an effect of dust obscuration 
by intervening galaxies [8]. 

All these hypotheses are indeed constrained by the recent observations of hight redshift quasars (z > 4) ; 
in fact, the effectiveness of the Gunn Peterson test increases with the redshift, the optical depth tgp being 
proportional to (1 + z)'^'^. 

Apparently, galaxies and unseen quasars cannot be the only photoionizing sources: in fact, the average 
ionizing intensity per unit of frequency and steradiant they produce is smaller that the value suggested by the 
so-called proximity effect; this effect is seen [9], [10] as a decrement in the counted number of absorbing Lya 
clouds in the immediate proximity of the known quasars: the Ly^ forest pattern is reduced as a consequence 
of the higher average ionization characterizing these clouds. 

The lower limit for the average intensity J_2i [9], [10] at the hydrogen Lyman limit (A = 912 A) is 
J_2i ^ 1 and it is independent on the redshift in the range 1.8 < 2; < 3.5; however, this lower limit for J-21 
is yet higher than the maximum flux coming from the known, counted quasars [10]. 

Thus, probably we need some additional radiation sources: in the following, I would like to discuss a model 
in which the evaporation of primordial black holes produces the high energy photons flux which ionizes the 
IGM at small redshifts. 

This model has been investigated in other papers: in particular, I discussed the effect of an exponential, 
late and fast reionization of the Universe on the polarization of the Cosmic Microwave Background in [11] 
while 1 discussed the details of the reionization processes induced by quantum evaporation of PBHs in [12], 
[13]. Here I will test the reliability of this model in the light of the Gunn Peterson effect, performing a 
comparison with the available experimental data. 

This paper is structured as follows: in Sect. 2 I discuss the main characteristics of the evaporating 
primordial black holes, particularly their mass evolution in presence of quarks and gluons jets emission; I 
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(1.1) 



also recall the main eqiiations giving the photon emission spectrum: more details about the time evolution 
of the ionization degree x and of the plasma temperature Tg can be found in [13] . 

In Sect. 3, I discuss the average ionization intensity J_2i coming from the photons emitted by PBHs during 

their evaporation: the lower limit J_2i ^ 1 suggested by the proximity effect is employed to estimate the 
average PBH relics mass, Mrei, left after the evaporation of these objects. The agreement between theory 
and experiment is only possible for a PBH average relics mass Mrei ~ 10~^^ g: that means we need a 
complete evaporation of these primordial objects, the most effective reionization being produced at the end 
of the PBHs life. A deeper investigation about this point is under study: in fact, the quantum gravity effects 
might change the BH evolution when its mass lowers under the Planck mass. 

In Sect. 3 I also present the main results of this work: I calculate the HI Gunn Peterson optical depth 
TOP, both in the case of a homogeneous IGM and in the more realistic case of an IGM presenting some 
inhomogeneities; as in ref. [8] , a clumping factor / takes into account the presence of moderate overdensities 
in the IGM. The cases of low (LA), medium (MA) and high (HA) photoelectric absorption by Lya clouds 
are all examined for different values of the density parameter 0/gm, namely fi/cM = 0.010, 0.015, 0.020: 
larger values produce a too high tgp that disagrees with the experimental data concerning high redshift 
{z ^ 4) quasars. 

In Sect. 4, 1 discuss in some detail the clumped structures present in the IGM in the form of Lya clouds, 
I determine the efficiency of ionization Gh and I present my results for the hydrogen number density nH,c, 
the related column density Nhi and the density parameter ClLa- 

Two important quantities are involved in this calculation: the Doppler shift parameter b and the cloud 
temperature T^, some evidence of a positive correlation b — Tc and b — Nhi was suggested in ref. [14] and 
it has been related to the presence of expanding, adiabatically cooled clouds in [15]; my analysis seems to 
confirm such a possibility. 

In Sect. 5, I perform the same calculation in the case of ionized helium and I discuss the relevance of 
the Gimn-Peterson effect for He II; the intergalactic helium should be present in the form of singly ionized 
He II rather than in the neutral form, proving again the high ionization status of the IGM at high redshift. 
Jakobsen et al. [16] recently claimed the observation of the He II 304 A Gunn-Peterson absorption from 
diffuse IGM at z ^ 3.2, with a total optical depth TGp(Hell) > 1.7; the theoretical value I calculated is 
smaller than this preliminary limit, namely equal to 0.76. 

Finally, in Sect. 6 I summarize the results obtained by using this model of PBHs-induced photoionization 
of the IGM and I present my conclusions. 



2. THE QUANTUM EVAPORATION OF 
PRIMORDIAL BLACK HOLES. 



The primordial black holes are very interesting structures, from a theoretical point of view; however, the 
processes that caused their formation and, in general, their overall evolution are poorly known; that despite 
of the relevant effects their quantum evaporation may have on the present status of the Universe. 

Following a widely accepted idea, the black holes formation should characterize the early instants of the 
Universe after the Big Bang; however, the nature of the phenomena acting to create such structures is not 
clear and many theories have been advanced [17]. 

The creation of a black hole is induced by the contraction of a mass to a size less than its gravitational 
radius; large mass stars at the end of their evolution may typically represent some good candidates to form 
BHs. Indeed, a gravitational contraction of such a size is quite unlikely for small stars and therefore the 
formation of black holes having a small mass (M < lO^'' g) is only possible in presence of a huge compression, 
i.e. at the early stages of the cosmological expansion. 

These small mass, primordial black holes present an intense Hawking's quantum evaporation: in fact, 
the blackbody temperature of the emission is inversely proportional to their mass [18] : 



kT = ——- ~ 1.06 

SttGM 



M 



-1 



10-, J ^'-'^ 

The nature of the emitted particles clearly depends on the blackbody temperature: the mass loss makes the 
BH hotter and hotter, thus enabling the emission of more and more massive particles. BHs having a mass 
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larger than 10^^ q can only emit masslcss particles [19], [20], [21]; however, when the BH mass falls below 
10^4 g, hadrons emission becomes possible: for temperatures above the QCD confinement scale Aqcd, a 
relevant emission of quarks and gluons jets does start; thus, the resulting spectrum is no more a blackbody 

one [20], [21], [22]. 

As I discussed in [12], [13], the photons emitted by PBHs during their quantum evaporation may be 
one of the causes of the reionization of the Universe: indeed, they seem very efficient to produce a late 
and sudden (nearly exponential) rise of the ionization degree for a reionization redshift zr < 30, while for 
30 < zn < 60 they are still able to cause a partial reionization. 

Note that the plasma heating due to the photons / electrons interactions is not strong enough to induce 
a relevant distortion of the CBR spectrum: that is in agreement with the recent FIRAS upper limit on the 
comptonization parameter, yc < 2.5 x 10~^ [23]. 

The most significant particle emission should happen at the last stages of the life of a BH; thus I choose 
a reionization time which corresponds to the evaporation time for an object having a mass near the critical 
one (Mc ~ 4.4 x lO^** h~°-^ g is the mass of a PBH that survives till the present epoch). 

In the following, I will briefly recall the fundamental equations describing the quantum evaporation of 
a PBH: more details can be found in [12], [13] and references therein. 

The initial mass Mj of a PBH is connected to its lifetime by the following general formula [22] : 



te 



1.19 X 10^ 



G^Mf 

nc^ f{Mi) 



6.24 X 10-2^ /(Mi)-i Mf sec; 



(2.2) 



the function f{M) in eq. (2.2) contains the contributions of the different species of particles, it is normalized 
to the unit for very massive (M > 10^'' g) BHs and reads as follows: 



f{M) 
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(2.3) 



In eq. (2.3) the first addendum in the right-hand side expresses the contribution of electrons, positrons, 

photons and neutrinos; heavier particles arc considered in the remaining terms, following their relative 
importance and with a factor 3 taking into account the color charge for the quarks; the denominators in the 
exponential terms are defined as the product PgjMj, where Mj is the mass of a black hole whose temperature 
is equal to the rest mass fXj of the j specie and Psj is a spin-dependent factor defined [22] in such a way the 
energy of a BH having M = (3sjMj has a peak at /Xj. 

When the BH mass falls below IQi^ g and the temperature T becomes larger than the confinement scale 
Aqcd, it is no longer possible to neglect the quarks and gluons emission in the calculation of the BH mass 
evolution: following [20], [21], [22], one should write this evolution as 



dM _ 1 



2tt% 







r SttGQM 1 


h 


exp 





Q dQ 



(2.4) 



in eq. (2.4), Tj is the absorption probability for the j particle having a spin Sj [24] and one sums on all the 
emitted species [22]; this equation means the emission of a parent particle j with total energy Q decreases 
the BH mass by Q/c^. After the integration over the energy Q, eq. (2.4) can be rewritten as [22]: 



dM 
dt 



— = -5.34 X 10^^ /(M) M 



^ g sec ^. 



(2.5) 
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Now, the Hawking emission rate of particles having an energy in the range (E, E + dE) from, a black hole 
having an angular velocity ui, an electric potential (j) and a surface gravity n is [18]: 



dN T dE E-nlt^-e<p\ , ^ 

where the signs ± respectively refer to fermions and bosons and F is the absorption probability of the emitted 
species. In the case oi photon emission, F reads [24]: 

- 1 i^r & <") 

in eq. (2.7) A is the surface area of the BH and the Planck mass and energy assure we are working with 
dimensionless quantities, as in [24]. 

Here, I neglect the charge and the angular momentum of PBHs, a quite reasonable and simplifying 
assumption due to the fact their loss via the quantum evaporation happens on a time scale shorter than the 
one characterizing the mass evaporation [25]. 

As pointed out in [26] in the emission spectrum we need a fragmentation function in order to take into 
account the production of quarks and gluons jets: 

dN^ V- /•+°° r,(Q,T) / Q ^,\-^ dgj^{Q,E) 



dtdE 



here x and j respectively label the final and the directly emitted particles while the last factor, containing 
the fragmentation function gj^, expresses the number of particles with energy in the range {E,E + dE) 
coming from a jet having an energy equal to Q [26]: 

in eq. (2.9) nih, is the hadron mass, fc is a constant ~ 0(1) and m is an index equal to 1 for mesons and 2 
for baryons. 

After selecting the dominant contribution to the integral over Q and summing over the final states, eq. 
(2.8) becomes [26]: 

^ E^ earp(^) for E»T Q E, (2.10a) 

E-^ far T~E»mh Q^T, (2.106) 



dtdE 
dN dg 



for E r.mh«T Q ~ m^; (2.10c) 



dtdE dE 

eqs. (2.10a), (2.10b) and (2.10c) respectively hold for the dominant value of Q written on the right. 



3 THE IONIZING PHOTON FLUX 



The possible origin of the ionizing flux has been studied in many papers [8], [27], [28], [29], [30], [31]. 
Generally, quasars are coinsidered the best candidate sources for the photoionization of the Universe; however, 
a background generated by quasars only is not consistent with the predictions of the proximity effect and 
with the data concerning the evolution of Lya clouds. For instance, in ref. [31] the Authors claimed that 
these quasars should ionize the IGM too late and produce a too large Gunn Peterson optical depth (note 
however that their conclusions have often been questioned). 
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In the following, I will discuss the mean specific intensity of the radiation field coming from the evapo- 
ration of PBHs as a function of the final mass of their relics; then, I will calculate the Gunn Peterson optical 
depth associated with the resonant HI Lya absorption. 



3.1 THE RADIATION INTENSITY 
AT THE HYDROGEN LYMAN EDGE 



At the hydrogen Lyman edge, ul = c/912 A = 3.29 x 10^^ sec ^, and for an observer at a redshift 
Zobs, the mean intensity of the photons flux in erg cm~^ sec~^ Hz~^ sr~^ is [32]: 



•^912(^065) = 



c (1 + Zob,f eriv^z) exp [-T^ff {912, Zp^s, z)] , ^ 



in eq. (3.1.1) erii^, z) is the total emissivity of a primordial population of BHs, calculated at a frequency 
u = (1 + -z) (1 + Zohs)^^ and at a redshift z, expressed in erg cm~^ sec~^ Hz~^; Teff {912, Zobs, z) is an 
eff'ective optical depth that, as I will discuss below, takes into account the probability that these photons are 
absorbed. 

The emissivity e can be obtained by calculating the total photon number density from eq. (2.10c). In 
our case, the condition E «T holds: we are effectively interested in the last stages of the evaporation, when 
the BHs temperature is very high; moreover, the processes that are relevant for this study (the ionization 
and the recombination) dominate [20], [21] for an energy E < lA KeV. 

For a jet fragmentation function given by eq. (2.9), the photon density in the unit time and energy can 
be written as follows: 



^ ^Tot = a. ZA mcs + ^jsor, (3.1.2) 



ocoot aivot ouiot 

where the mesons and baryons contributions have been singled out. 
For Q = ruhadr ~ 300 MeV, eq. (3.1.2) reads 

the dominating contribution being the mesonic one; after reducing to the proper dimensions, I obtain the 
volume emissivity of a single evaporating PBH: 

dn^PR , X 3.20 X 10-15 3 i ,-r i /o -■ .x 

OTOw w (Ge v ) 

After considering the frequency shift y = vl (1 + -J) (1 + Zobs)~^ , eq. (3.1.4) becomes: 

^ 3-^6 X 10-^0 {1 + Zobs) ^^_3 (3_^_5^ 

^ ' VL {1 + z) ^ ^ ' 

Finally, the total emissivity is: 

eT{i',z) = Npbh{z) e{iy,z), (3.1.6) 

where the parameter Npbh{z) represents the population of PBHs at a redshift z. 
Its value can roughly be estimated as follows: I write the PBHs density as 



Mr,l{z) Npbh{z) 

"^'^ W{zj ' ^^-^-^^ 

where Mrei{z) is the average mass of the PBHs relics at a redshift z. Then, the scale factor evolves as in 
the radiation dominated epoch: 

R{t) ~ Ro {t/tof^; (3.1.8) 



{Rq = 1.25 X 10^** cm = 1.4 x 10^" lyr in a standard cosmological model, sec rcf. [33]). 

Now, I assume the PBHs density evolution is approximately described by a power law with index 2/3; 
the formation time coincides with the Big Bang and the initial density is p, ^ 4.28 x 10^^ g cm~^, [13] as one 
infers on the basis of the present experimental limits on flpBH {'^pbh < (7.6 ± 2.6) x 10~® ^(-i.95±o.i5)'j 
[20], [21]. Then, one finally gets [13]: 

Npbh{z) = Ro Pi C^f (1 + z)'/' , (3.1.9) 

to M rel [Z ) 

i.e.: 

NpBHiz) = 9.67 X 10-^'' ^-'^^^ cm-2. (3.1.10) 

Mrel{z) 

The most difficult parameter to estimate in eq. (3.1.10) is the average relics mass Mrei at the reionization 
epoch {z ~ 20-f30); 

This parameter can be evaluated by searching a basic agreement with the data coming from the proximity 
effect and concerning the average intensity of the ionizing background; an agreement is only possible for 
Mrei ~ 10~^* g, i.e. for black holes totally evaporated or may be for relics well lighter than the Planck mass. 
However, for these objects the quantum gravity effects might be very relevant, thus a deeper theoretical 
analysis about this point is under study. 



Turning back to the effective optical depth, it can be written as [34]: 



Teff{912,Zobs,z) = 



Zabs JO 



dNuidz 



7 X [1 - exp{-NHi(Tv')] dNni dz ; 



(3.1.11) 



in eq. (3.1.11) a^,' is the hydrogen photoionization cross section, Nhi is the hydrogen column density of the 
absorber and the derivative dN/dz represents the column density distribution, depending on the assumed 

model of attenuation [32]: in fact, we can have various level of absorption and the approximate integration 
of eq. (3.1.11), respectively in the cases of low, medium and high absorption, gives [32]: 



T^f^f{912,Zobs,z) ~ [ 0.0118 a;^,^ {x''-^ - x^obs) + 2-35 ar^.^, ln( — )- 

-0.78 {x-Js"" - X-'-') - 0.003 {x'-' - xl,l) ] , 

Ti^/(912,^„,„^) ^ \0.2Ux%Ax'-'-xl,t) + 2.35 xliln( — 
-0.78 {x-Jf - x-i-S) - 0.003 (xi-5 - x]^l) 



T^fj [912, Zobs,z) 



0.097 xii,^ (x^ «^ - xi(,r) - 0.0068 x^,, {x 



0.4 



Xobs) 



-8.06x10-^ {x^-^-xl,\) 



(3.1.12) 



(3.1.13) 



(3.1.14) 



where Xobs = 1 + Zobs and x = 1 + z. 

The effective nature of the absorbers is still unclear: probably, they are large clouds (up to 1 Mpc) 
containing a significant baryonic fraction [35]; many recent analyses identify such absorbers with evolved 
densities fluctuations in the intergalactic medium, produced by large scale flows and inhomogeneities of the 
dark matter component of the Universe [36] ; anyway, these absorbers both severely attenuate the flux coming 
from the ionization sources and generally delay the growth of HII regions [32]. 

By using eqs. (3.1.5), (3.1.6) and (3.1.10) in eq. (3.1.1), I obtain: 

C (l + ZobsV 



J912{Zobs) 



pzmo 
J zobs 



(1 + 2)25/4 



exp [-Tef f {912, Zobs, z)]- 
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9.67x10-1° 3.36x10"*° , 
• = dz. (3.1.15) 

Mrel 

After recalling the definition of the reduced intensity 

J-2l{Zohs) = JMZobs)/10~^\ (3.1.16) 

I plot in fig. 1 the behaviour of J-2i(zobs) vs Zobs obtained from eq. (3.1.15), in the cases of low (LA), 
medium (MA) and high (HA) absorption; I take Mrei = Mrei{zreion) ~ 10""'^^ g. 

The Gunn-Peterson optical depth associated with the resonant Lya absorption can be written in the 
following form [1]: 

^ \ - ( ^ nHi,d{z) ^'^^^7\ 

TGP[Zobs) - I TT ItT- , „ 7T72> (3.1.17) 

where nHi,d is the HI density of the IGM diffuse component and fa is the oscillation strength of the Lya 
transition. 

Rewriting eq. (3.1.17) as a function of the intensity J-2i{zobs) one obtains [32]: 

TGp{Zobs) = 3 T4-0-^5 (^^^^ /^2)2 (3 ^ {1 + z)^'' f{z) {J -2l{Zobs))~^ \ (3.1.18) 

eq. (3.1.18) is evaluated by inserting the temperature T4 = T/10^ K° of the ionized intercloud medium, as 
it results from the PBHs induced photoionization processes [12], [13]. A plot of T4 (in K°) as a function of 
the rcdshift z is shown in fig. 2; in eq. (3.1.18) a is the power spectrum index of the metagalactic flux at 
high z {J (X in our case 

J-2l{Zobs,v) = -^-21 (^obs), (3.1.19) 

thus a = 1. 

If compared to the original equation of Gunn and Peterson, Eq. (3.1.18) contains an additional factor, 
namely the clumping factor f{z): the presence of inhomogcncitics in the IGM may influence in a relevant 
way the processes of absorption [5], [6], [7], [8], [37] and, indeed, the observation of Lya clouds is a clear 
proof of the existence of remarkable over densities. 

The function / describes moderate overdensities (1 < p/p < 10) and it is deflned as follows [8]: 

10 - /, (3.1.20) 

while larger inhomogcncitics arc directly identified with Lya clouds. 

The values of / at some different redshifts are listed in tab. 1; these values have been obtained by 
Carlberg and Couchman [38] through a numerical simulation. 

Another important parameter entering in eq. (3.1.18) is the density parameter for the intcrgalactic 
medium, ^igm'- here I tested 3 values, namely ^igm = 0.010, 0.015, 0.020, for all the cases of low, medium 
and high absorption and both for an idealized, perfectly homogeneous IGM {f{z) = 1) and in presence of 
moderate inhomogcncitics {f{z) as in tab. 1). Figs. 3a, 3b, 4a, 4b, 5a, 5b show the behaviour of tgp vs 
Zohs for all these cases; in tab. 2a, 2b, 3a, 3b, 4a, 4b I listed the values of tgp calculated for two reference 
redshifts, Zobs = 3, z^bs = 4-3, for which some experimental data are known; tab. 5 finally gives a picture of 
the present experimental tqp measurements, that also includes the recent data of Giallongo ct al. [39]. 

Finally, in figs. 6a, 6b and 6c I resume the behaviour of tgp{z) for the different values of the density 
parameter fljGM and for various levels of absorption. 

3.2 DISCUSSION OF THE RESULTS 

As I told, the average photon intensity J-21 determined from the proximity effect is employed to 
constrain the final average mass of the PBHs relics, that should be nearly totally evaporated at the time of 
the reionization. 
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The calculation of the Gunn Peterson optical depth is strongly model-dependent: in fact, eq. (3.1.18) 
contains both the temperature of the ionized IGM and the power spectrum index a of the metagalactic flux 
at high z, i.e. two peculiar predictions of the model. 

Looking firstly at fig. 1, one can remark that the radiation intensity evolution depends on the behaviour 
of the effective optical depth, eq. (3.1.11); the same behaviour is obtained for J_2i in the cases of low and 
medium absorption, only rescaled by a factor weighting the first cubic term in eqs. (3.1.12) and (3.1.13). 
On the contrary, in the case of a high absorption level, eq. (3.1.14) for r^ff produces a steeper curve. 

The behaviour of J_2i vs. z differs in a sensitive way from the one predicted, for instance, in [8], 
where the sources of the ionization one assumes are quasars and high-mass stars in primordial galaxies. At 
early epochs, the intergalactic gas is neutral and the Universe is opaque to the radiation, while the volume 
emissivity of photons, s{i', z), is not so high to produce a relevant ionization; thus, the intensity J_2i in fig. 
2 goes to zero at high redshifts. 

In ref. [8] the rise of the intensity J_2i begins at z ~ 5 and stops at z ^ 2.5, due to fact that, at this 
redshift, the quasars start to decline; on the contrary, the behaviour of the intensity shown in fig. 1 is nearly 
exponential, due to the peculiar emission process I considered. 

Looking at tabs. 2a, 2b, 3a, 3b, 4a and 4b, together with the experimental available data for the Gunn 
Peterson optical depth shown in tab. 5, we can observe that: 

a) the experimental data generally refer to different values of the redshift z and therefore, in order to allow 
a comparison, I considered two reference redshifts z = 3 and z = 4.3 for which many measurements 
have been performed. 

b) The clumping factor / enters in a linear way in eq. (3.1.18): thus, the presence of some inhomogeneities 
(/ > 1) in the intergalactic medium increases the GP optical depth; a clumped configuration of the 
IGM makes it more opaque to the radiation. 

c) I tested three values of the IGM density parameter, i.e. flicM = 0.010, 0.015, 0.020; values larger 
than 0.020 produce a too high optical depth top, disagreeing with the results of Giallongo et al. [39] 
(that however are affected by significant uncertainties). The disagreement is particularly serious in the 
inhomogeneous case (/ > 1) with a medium / high absorption. 

I can obtain a satisfactory agreement with the data of Giallongo et al. [39] for ilioM = 0.010 and a low 
absorption level: smaller values of fi/cM are in principle acceptable but, as I will show, they do not allow 
to obtain a consistent value of the diffuse medium density parameter, Qd = ^igm — ^Lya- 

4 CHARACTERISTICS OF THE Ly^ CLOUDS 

4.1 THE THEORY 

Here I want to discuss the characteristics of the IGM regions presenting very large overdensities, p/p > 
10, i.e. the configurations known as Lya clouds; many informations can be obtained by studying the Lya 
absorption phenomena within these regions. A very interesting idea is the possibility that these clouds are 
expanding and adiabatically cooled: the observations of Pettini et al. [14] suggested that the Ly^ forest lines 
with a low column density Nhi {Nhi < 10^^ cm~^) have also very small velocity widths b. They found a 
positive correlation between b and Nhi; typical values for these parameters are [14], [15]: 

< 6 > ~ 11 ± 3 /cm sec-^ for log Nhi = 13; (4.1.1) 

Nhi is expressed in cm~^ and b (also called " Doppler parameter") is expressed as 

= Vl,, + ^, (4.1.2) 

i.e. it is connected to the cloud characteristics, as the temperature and the bulk velocity. 

Generally, these clouds are considered as small, dense and very filamentous structures, rather than 
large and highly ionized objects; however, a different interpretation has been proposed in [15], [40] where 
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one assumes that these low density clouds are indeed expanding structures, which progressively cool in 
an adiabatic way; the expansion cooling may bring highly ionized clouds out of the thermal ionization 
equilibrium, thus enabling lower values of temperature (T << 10^ K°) and of the velocity widths but not 
necessarily a low ionization level. 

In the following, I will investigate the characteristics of Ly^ clouds by studying their absorption of the 
PBH photons; in particular, I will calculate the ionization efficiency, the HI column density and the baryonic 
content of the Lj/q clouds, expressed by the density parameter ^Lya- We will sec that, in such a model 
of PBH-induced reionization, senseful results for Q. can be obtained only in presence of expanding and 
adiabatically cooled clouds; thus, my results seem to support the hypotheses of ref. [15]. 



4.2 IONIZATION EFFICIENCY, HI DENSITY 
AND DENSITY PARAMETER QLya 



Here I will recall the main equations that are useful to study the ionization problem. In the case of 
photoionization equilibrium between the cloud gas and the ionizing radiation and for a ionizing flux given 
by eq. (3.1.15), the ionization efliciency is [32]: 

Gh = r ^ ^ ^, (4.2.1) 



hp V J912 



where hp is the Planck constant, JvlJ^n is given by eq. (3.1.19) and the hydrogen photoionization cross 
section is [41] 



(Jy = An 



( VL\^ exp [4- ((4arctan£)/£)] 
°IVJ l-exv (-27r/e) ' ^^-^"^^ 



holding for v = vl and with 



1 — exp {—2-k/e 



^0 = {^) -^rl = 6.30 X 10-i« cm^ (4.2.3) 



— - 1. (4.2.4) 



Eq. (4.2.1) can be recast in the simpler form 



Gh ~ 2.11 X 10" UL I ^ du, (4.2.5) 



where a suitable upper cut in the integration {vmax = 1. x 10^^ sec ^) is inserted in order to numerically 
handle the integration, after testing that this choice does not seriously influence the final result. I obtain: 

Gh = 3.52 X lO'^^ erg''^ sec. (4.2.6) 

The total hydrogen gas density Uh^c of a spherical cloud of average column density Nhi = 10^'' cm~^ is 
related to its diameter D in the following way [32]: 



nH,c{z) = 



3 GgJ9i2(z)(10^'^ cm-^ 
2 (l + 2x) aA{Tc) D 



-,1/2 

(4.2.7) 



where aA = 4.2 x 10 ^'^ T4 cm? sec ^ is the coefficient of recombination to all the levels of hydrogen 
[41], T4_c is the cloud temperature {K) divided by 10** and x is the ratio He/H, equal to 1/12. The column 
density Nhi, averaged over all the lines of sight, is connected to the diameter D of the cloud through the 
following formula [32] : 

<Nhi>= lnH,cD, (4.2.8) 
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and the Lya clouds density parameter, O.Lya is 



^Lya = fc (4.2.9) 
nH,crit 



here the closure hydrogen density is 



and the volume filling factor /c of the Ly^ forest is [32]: 

/e(^) ~ D{z) (1 + z)^-^ (1 + 290^)°-^ (4.2.11) 
Finally, the size of the Lya clouds is connected to the Doppler parameter b by the approximate formula [32]: 

D{z) ~ < 6 > [Ho{l + ^) (1 + 2qozf^]-\ (4.2.12) 

Choosing the standard values Tc = 2 x 10* K° and < b >= 35 km sec~^ [15], I obtained for the cloud 
hydrogen gas density uh.c and for the column density Nhi the values listed in tab. 6a for two reference 
redshifts, z = 2.5 and z = 0.0. 

The baryonic density parameter at the present time for the Lya clouds I calculated from eq. (4.2.9) is 

QLya = 0.016; (4.2.13) 

This high value means that a large fraction of baryons resides in the Ly^ clouds, while the diffuse component 
of the IGM = ^IGM — ^Lya) should have a density parameter totally negligible. A more reliable 

result can be obtained if we consider the clouds as dynamically expanding structures, a possibility that also 
could explain the (probably) observed correlation between the Doppler parameter b and the column density 
Nhi. 

In this case, the temperature Tc lowers to 200 K° and the size grows up to 1 Mpc; the results I obtain 
are listed in tab. 6b. For this choice of Tc and D, the present density parameter is 

^Lya = 0.010, (4.2.14) 

in all the cases of low, medium and high absorption (see tab. 7). Consequently, the diffuse medium should 
have a present density parameter 

< fir, < 0.010, (4.2.15) 

depending on the effective value chosen for flioM- 

Above, I calculated the hydrogen density for the clumped component of the IGM; now, the Gunn 
Peterson optical depth tgp{z) theoretically obtained in the previous section, can be used in eq. (3.1.17) in 
order to predict the density nni. d of its diffuse component. 

Eq. (3.1.17) can be rewritten as follows: 

8.28 X 101° UHi, d{z) 

= OTIF^ ' ^'-'-'^^ 

and therefore: 

UHi, d{z) = 1.21 X 10-" TGp(z) (1 + zf/\ (4.2.17) 

The density nni, d calculated for three reference redshifts, ^ = 0, z = 2 and z = 2.64 are listed in tabs. 8a, 
8b and 8c for CIigm = 0.010, 0.015, 0.020 and for the cases of low, medium and high absorption. 
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5. THE HE II GUNN PETERSON TEST 



An interesting extension of the Gunn Peterson test is the study of the Hell 304 A resonance hne of the 
singly ionized helium; both the theoretical and the experimental efforts are converging on this problem since 
the appearance of the results of Jakobsen et al. [16], concerning the probable detection of the Hell 304 A 
Gunn Peterson absorption from the IGM at ^; 3.2. 

The primordial helium in the intergalactic medium should be strongly ionized at high redshifts and it 
should be mainly seen in the form of singly ionized Hell, rather than neutral Hel [16]; this fact has been 
confirmed by the failure in detecting the Gunn Peterson absorption in the Hel 584 A line of the neutral 
helium [42]. 

The Hell Ly^ resonance line is redshifted at Xobs ~ 304 {1 + z) A and therefore it is observable in the 
extreme UV; thus, the He II Gunn Peterson test could potentially give a lot of informations on the spectrum 
of the UV background at very short wavelengths. Distant quasars at redshifts z > 3 are the most useful 
objects to study, in order to detect this effect in an experimentally suitable range of wavelengths, A > 1200 A, 
where the absorption by hydrogen in our galaxy does not produce noise. 

Prom a theoretical point of view, the same kind of analysis used for the hydrogen can be performed for 
Hell, with only few differences; firstly, the Hell ionization edge is located at A = 228 A, corresponding to a 
frequency vl, He = 1-31 x 10^^ sec; from eqs. (3.1.17) and (3.1.18) it is possible to estimate [32] the ratios 
of the intensities of the diffuse radiation field at the hydrogen Lyman edge {i/l, h = c/912 A) and at the 
helium edge [ul, He = c/228 A); for a low attenuation model, Meiksin and Madau found [32]: 

-^912 / -^228 ~ 8, (5.1) 

while in a medium attenuation model [28] one approximately has 

■^12 / -^228 ~ 25. (5.2) 

From eqs. (5.1) and (5.2) I simply obtained the intensity J22&, that I plotted in fig. 7. 

Using those intensities in eq. (3.1.18) I obtain the results shown in tabs. 9a, 9b, 10a, 10b and 
11a, lib for the Gunn Peterson optical depth at two reference redshifts, z = 3 and z = 4.3 and for 
^iGM = 0.010, 0.015, 0.020. A plot of top{z) for the values of flicM and both in the homogeneous 
and inhomogencous cases is shown in figs. 8a, 8b. 9a, 9b, 10a and 10b; finally, in figs. 11a and lib I resumed 
the behaviour of tgp{z) for all the values of ilicM and for a fixed attenuation level. 

The analysis of §4.2 can be repeated, simply observing that the Lya clouds absorption process is just 
modified as follows: the He absorption cross section is connected to the H cross section through the relation 
[8]: ^ 

(THeIl{v) = ^ CTff (^); (5-3) 

The ionization efficiency Gne can be simply calculated by using eq. (5.3) in eq. (4.2.1) and changing the 
intensities: thus, I calculate 

GHe = 2.25 X 10"^° erg-'^ sec, (5.4) 

Finally, the helium density uhc, c and the column density Nneii are given by some equations analogous to 
eqs. (4.2.7), (4.2.8), with a ionization efficiency given by eq. (5.4), an intensity J228 coming from eqs. (5.1) 
(LA case) and (5.2) (MA case) and a coefficient of recombination for He given by [41]: 

aA{T) = 2.607 x 10"" T^° ''^. (5.5) 

The results are listed in tab. 12a for Tc = 2 x 10^ K° and in tab. 12b for adiabatically expanding, cooled 
clouds at Tc = 200 K°. 

A comparison of the Gunn Peterson optical depth here obtained with the experimental value of Jakobsen 
et al. [16] {THeii > 1-7 at 2; ~ 3) is also possible but this lower limit should be considered with some prudence: 
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indeed, it is very conservative and probably affected by some unc;c;rtainties, due to the fact that the Authors 
cannot distinguish between Hell absorption coming from the hne blanketing in the discrete lines of the Lya 
forest and the real Gunn Peterson trough. 

Looking at tabs. 9a, 9b, 10a, 10b and 11a, lib we can remark that in such a PBH-induced model of 
reionization, the Gunn Peterson optical depths obtained for Hell are quite similar to the ones for HI: in the 
case of inhomogeneous clouds with a medium absorption level and for CIigm = 0.020, the maximum value 
for Top hcII at z = 3 is 0.76. The disagreement with the lower limit of ref. [16] should be tested in the light 
of future, more precise measurements. 

6. CONCLUSIONS 

In this paper I assumed that the quantum evaporation of Primordial Black Holes having a mass M ~ 

10"'^^ g provides the radiation that photoionizes the intergalactic medium at some early epochs in the past. 

This particular model has been tested in the light of the Gunn Peterson test for both HI and Hell: I 
calculated the optical depths associated with the resonant Lpa absorption by clouds and I compared my 
theoretical resiilts with the last experimental data. The agreement is satisfactory, particularly if we assume 
an inhomogeneous IGM with a density parameter Q = 0.010 and a low attenuation by the clouds. In any 
case, a constraint on Q.igm can be inferred from this analysis: the agreement theory / experiment is possible 
only if niGM < 0.020. 

I also studied the characteristics of the Lya clouds and I found that my model confirms the possibility 
they are expanding and adiabatically cooled. 

A comparison theory / experiment for the Hell Gunn Peterson effect is certainly premature, due to the 
paucity of the experimental results regarding Hell: the Hell optical depth I calculated is smaller than the 
lower limit claimed by Jakobsen et al. [16] [tgp > 1-7) and near to the values I found for HI. Finally, I can 
conclude that the reionization model here proposed seems to work quite well; anyway, a better knowledge of 
the ultimate fate of the PBHs relics, particularly of their final masses, surely represents a definitive test for 
it. 
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Tab. 1: Clumping function f coming from a numerical simulation of Cctrlberg and Couch- 
man (Ccirlberg & Couchman 1989) 



f 


z 


7.42 


0.00 


7.98 


0.81 


6.93 


1.36 


2.48 


2.80 


1.74 


4.63 
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Tab. 2a: Gunn Peterson optical depth: homogeneous case (/ = 1) with a density parameter 

^IGM = 0.010 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 


0.002 


0.014 


Medium Abs. 


0.003 


0.030 


High Abs. 


0.003 


0.036 



Tab. 2b: Gunn Peterson optical depth: inhomogeneous case (/ > 1) with a density 
parEimeter 0/gm = 0.010 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 


0.005 


0.026 


Medium Abs. 


0.008 


0.057 


High Abs. 


0.008 


0.068 
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Tab. 3a: Gunn Peterson optical depth: homogeneous case 
(/ = 1) with a density parameter iljcM = 0.015 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 


0.004 


0.031 


Medium Abs. 


0.007 


0.068 


High Abs. 


0.007 


0.081 



Tab. 3b: Gunn Peterson optical depth: inhomogeneous case (/ > 1) with a density 
parameter 0/gm = 0.015 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 


0.010 


0.057 


Medium Abs. 


0.017 


0.128 


High Abs. 


0.017 


0.152 
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Tab. 4a: Gunn Peterson optical depth: homogeneous case 
(/ = 1) with a density parameter iljcM = 0.020 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 


0.008 


0.055 


Medium Abs. 


0.012 


0.121 


High Abs. 


0.013 


0.144 



Tab. 4b: Gunn Peterson optical depth: inhomogeneous case (/ > 1) with a density 
parameter 0/gm = 0.020 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 


0.018 


0.102 


Medium Abs. 


0.030 


0.230 


High Abs. 


0.030 


0.270 
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Tab. 5: Gunn Peterson optical depth: experimental data 



TGP 


z 


year 


Work 


0.01 


3.0 


1992 


Giallongo et al. " 


0.01 ±0.03 


3.0 


1994 


Giallongo et al. ^ 


0.02 ± 0.03 


4.3 


1994 


Giallongo et al. * 


< 0.05 


3.8 


1992 


Webb et al. " 


0.04 


4.0 


1992 


Webb et al. " 


0.04 ±0.01 


4.1 


1992 


Webb et al. " 


< 0.02 ± 0.03 


2.64 


1994 


Steidel and Sargent 



a) Giallongo et al., APJ 398, L12, 1992. 

b) Giallongo et al., APJ 425, LI, 1994. 

c) Webb ct al., MNRAS 255, 319, 1992. 

d) Steidel and Sargent, APJ, 318, Lll, 1987. 
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Tab. 6a: Lya clouds characteristics: hydrogen density Uhc scad column density Nhi at a 
reference redshift z = 2.5 and for Tc = 2 x 10^ K° 



ABS. MODEL 


hhc cm ^ 


Nhi cm ^ 


Low Abs. 
Medium Abs. 
High Abs. 


1.04 X 10-"^ 
8.56 X 10-5 
8.91 X 10-5 


2.28 X 10^9 
1.88 X 10^9 
1.96 X 10^9 



Tab. 6b: Ly^ clouds characteristics: hydrogen density uhc and column density Nhi at a 
reference redshift z = 2.5 and for Tc = 200 K° 



ABS. MODEL 


hhc cm ^ 


Nhi cm ^ 


Low Abs. 
Medium Abs. 

High Abs. 


1.84 X 10-5 
1.52 X 10-5 

1.58 X 10"''^ 


4.05 X 10^^ 
3.35 X 10^^ 

3.48 X IQi* 
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Tab. 7: Density parameter flLya for Lya clusters for a cluster temperature Tc = 200 K° 



ABS. MODEL 


= 2.5 


= 0.0 


Low Abs. 


0.013 


0.010 


Medium Abs. 


0.011 


0.010 


High Abs. 


0.012 


0.010 



21 



Tab. 8a: Hydrogen density nn. d at redshift z = 0, z = 2.00 and z = 2.64 characterizing the 
diffuse component of the IGM; ^igm = 0.010 and / > 1 



ABS. MODEL 


nH,d{z = 0) cm ^ 


nH,d{z = 2-00) cm-3 


nH,d{z = 2.64) cm-3 


Low Abs. 
Medium Abs. 
High Abs. 


4.91 X 10-1^ 

5.01 X 10-1^ 

4.02 X 10-1^ 


8.99 X 10"" 
1.19 X 10-13 
1.01 X 10-13 


2.47 X 10-" 
3.83 X 10-13 
3.62 X 10-13 
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Tab. 8b: Hydrogen density uh, d at redshift z = 0, z = 2.00 and z = 2.64 characterizing the 
diffuse component of the IGM; SIigm = 0.015 and / > 1 



ABS. MODEL 


nH,d{^ = 0) cm ^ 


nH,d{z = 2-00) cm-3 


"ij,d(-z = 2.64) cm-3 


Low Abs. 
Medium Abs. 
High Abs. 


1.10 X 10-1^ 
1.12 X 10-1*^ 
9.04 X 10-1^ 


2.02 X 10-13 
2.67 X 10-13 
2.28 X 10-13 


5.71 X 10-13 
8.62 X 10-13 
8.15 X 10-13 
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Tab. 8c: Hydrogen density uh, a at redshift z = Q, z = 2.00 and z = 2.64 characterizing the 
diffuse component of the IGM; ^igm = 0.020 and / > 1 



ABS. MODEL 


nH,d{z = 0) cm ^ 


nH,d{z = 2-00) cm-3 


nH,d{z = 2.64) cm-3 


Low Abs. 
Medium Abs. 
High Abs. 


1.96 X 10-1^ 
1.61 X 10-1*^ 
2.03 X lO-^'^ 


3.60 X 10-13 
4.06 X 10-13 
4.75 X 10-13 


1.01 X 10-12 
1.45 X 10-12 
1.53 X 10-12 
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Tab. 9a: Hell Gunn Peterson optical depth: homogeneous case (/ = 1) with a density 
parameter ^igm = 0.010 



ABSORPTION MODEL 


z = Z 


z = 4.3 


Low Abs. 
Medium Abs. 


0.037 
0.190 


0.204 
1.420 



Tab. 9b: Hell Gunn Peterson optical depth: inhomogeneous case (/ > 1) with a density 
parameter ^igm = 0.010 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 
Medium Abs. 


0.015 
0.079 


0.109 
0.758 
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Tab. 10a: Hell Gunn Peterson optical depth: homogeneous case (/ = 1) with a density 
parameter Q.igm = 0.015 



ABSORPTION MODEL 


z = Z 


z = 4.3 


Low Abs. 
Medium Abs. 


0.083 
0.428 


0.460 
3.195 



Tab. 10b: Hell Gunn Peterson optical depth: inhomogeneous case (/ > 1) with a density 
parameter ftjcM = 0.015 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 
Medium Abs. 


0.035 
0.178 


0.246 
1.706 
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Tab. 11a: Hell Gunn Peterson optical depth: homogeneous case (/ = 1) with a density 
parameter Q.igm = 0.020 



ABSORPTION MODEL 


z = Z 


z = 4.3 


Low Abs. 
Medium Abs. 


0.148 
0.760 


0.818 
5.679 



Tab. lib: Hell Gunn Peterson optical depth: inhomogeneous case (/ > 1) with a density 
parameter fljcM = 0.020 



ABSORPTION MODEL 


z = 3 


z = 4.3 


Low Abs. 
Medium Abs. 


0.062 
0.317 


0.437 
3.032 
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Tab. 12a: Lya clouds characteristics: helium density URec and column density Nneii at 
reference redshift z = 2.5 and for Tc = 2 x 10'' K° 



ABSORPTION MODEL 


uhc cm ^ 


Nhi cm ^ 


Low Abs. 
Medium Abs. 


4.11 X 10-5 
1.92 X 10-5 


9.06 X lO^s 
4.23 X 10^^ 



Tab. 12b: Lya clouds characteristics: helium density UHec and column density NhbII at 
reference redshift z = 2.5 and for Tc = 200 K° 



ABSORPTION MODEL 


u hc cm. 


Nhi cm, ^ 


Low Abs. 
Medium Abs. 


2.32 X 10-'' 
1.08 X 10-^ 


5.09 X lO^^ 
2.38 X 10^9 
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FIGURE CAPTIONS: 



Fig. 1: Average ionizing intensity per unit of frequency and steradiant Jg2i{z) at the hydrogen Lyman edge 
ul = c/912 A {erg cm~^ sec~^ Hz~^ sr~^; full line: low attenuation model; dashed line: medium 
attenuation model; dotted line: high attenuation model. 

Fig. 2: Evolution of the temperature of the intergalactic medium as a function of the redshift z (eV). 

Fig. 3a: Gunn Peterson optical depth Tap{z) associated with the resonant HI Lya absorption for an inhomo- 
geneous IGM with Qigm = 0.010 and with a clumping factor / < 1; full line: low attenuation model; 
dashed line: medium attenuation model; dotted line: high attenuation model. 

Fig. 3b: Gunn Peterson optical depth tgp{z) associated with the resonant HI Lya absorption for a homogeneous 
IGM with fliGM = 0.010 and with a clumping factor / = 1. 

Figs. 4a-4b: The same plots of figs. 3a and 3b but for fljGM = 0.015. 

Figs. 5a-5b: The same plots of figs. 3a and 3b but for ^igm = 0.020. 

Fig. 6a: Gunn Peterson optical depth associated with the resonant HI Lya absorption for the three values of 
SljGM in the case of an inhomogeneous IGM (/ > 1) with low absorption; full line: flicM = 0.010, 
dashed line: f^/cM = 0.015, dotted lino: O/cAf = 0.020. 

Fig. 6b: The same plot of fig. 6a but in the case of medium attenuation. 

Fig. 6c: The same plot of fig. 6a but in the case of high attenuation. 

Fig. 7: Average ionizing intensity per unit of frequency and steradiant J228(-z) at the Helium Lyman edge 
i/L = c/228 A {erg cm~^ sec~^ Hz~^ sr~^; full line: low attenuation model; dashed line: medium 
attenuation model; dotted line: high attenuation model. 

Fig. 8a: Gunn Peterson optical depth tgp{z) associated with the resonant Hell Lya absorption for an inhomo- 
geneous IGM with iliGM = 0.010 and with a clumping factor / < 1; full line: low attenuation model; 
dashed line: medium attenuation model. 

Fig. 8b: Gunn Peterson optical depth tgp{z) associated with the resonant Hell Ly^ absorption for a homoge- 
neous IGM with 0/GM = 0.010 and with a clumping factor / = 1. 

Figs. 9a-9b: The same plots of figs. 8a and 8b but for CIigm = 0.015. 

Figs. lOa-lOb: The same plots of figs. 8a and 8b but for Qigm = 0.020. 

Fig. 11a: Gunn Peterson optical depth associated with the resonant Hell Lya absorption for the three chosen 
values of Hjgm in the case of an inhomogeneous IGM (/ > 1) with low absorption; full line: fliGM = 
0.010, dashed line: fljGM = 0.015, dotted line: fljGM = 0.020. 

Fig. lib: The same plot of fig. 6a but in the case of medium attenuation. 
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